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The ab-plane NIS-tunnelling conductance in d-wave superconductors shows a zero-bias conduc- 
tance peak which is predicted to spht in a magnetic field. In a pure d-wave superconductor the 
splitting is linear for fields small on the scale of the thermodynamic critical field. The field depen- 
dence is shown to be nonlinear, even at low fields, in the vicinity of a surface phase transition into 
a local time-reversal symmetry breaking state. The field evolution of the conductance is sensitive 
to temperature, doping, and the symmetry of the sub-dominant pairing channel. 
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Tunnelling into an unconventional superconductor 
probes the quasiparticle spectrum associated with a sur- 
face superconducting state that is typically strongly de- 
formed, and may have different local symmetry, than the 
bulk order parameter (OP) P|. Indeed, the OP of a 
pure dx2_y2 superconductor is suppressed at a [110] sur- 
face. The suppression is associated with the formation 
of zero-energy Andreev states that are bound to the sur- 
face (ABS). The ABS is observable as a zero-bias peak 
(ZBCP) in the tunnelling conductance This prop- 

erty makes quasiparticle tunnelling in unconventional su- 
perconductors a sensitive probe of broken symmetry. A 
key signature of a surface ABS is its dependence on an ap- 
plied magnetic field, or more precisely on screening cur- 
rents flowing along the surface. For an external magnetic 
field perpendicular to the conducting plane and parallel 
to the surface, screening currents produce a Doppler shift 
the ABS which is observable as a splitting, 25{H), of 
the ZBCP with applied field. Some of the predicted field 
dependence of the splitting has been verified experimen- 
tally by Covington et al. p, Aprih et al.Q, and Krupke 
and Deutscher for YBaaCugOr-p (YBCO) near opti- 
mal doping. 

It was shown theoretically that a surface phase transi- 
tion to a state with spontaneously broken time-reversal 
(T) symmetry should occur for a d-wave superconductor 
with strong surface pair breaking if there is an attractive 
subdominant pairing channel 0, I3- A signature of 
this transition in the tunnelling conductance is a sponta- 
neous splitting of the ZBCP that develops below a surface 
phase transition temperature, Tg . Such a transition was 
observed by Covington et al. l5j for Cu 1 1 1 YBCO junctions, 
and also by Krupke, et al. and Dagan et al. for 
In 1 1 1 YBCO junctions. The latter authors found that the 
spontaneous splitting of the ZBCP occurred only above 
a critical doping level, Pc, close to optimal doping. They 
also showed that the magnitude of the field- induced split- 
ting of the ZBCP is suppressed at low fields for slightly 
underdoped materials [p < pc), while for p w pc 6{H) 
is nonlinear in H and onsets rapidly. For slightly over- 



doped YBCO (p > Pc) S{H) increases linearly with H at 
low field from a non-zero splitting in zero field. 

In this Letter we report calculations of the field de- 
pendence of the conductance of NIS tunnel junctions in 
YBCO. We consider a pairing interaction that includes 
both mag netic and phonon-mediated pairing channels 
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The repulsive magnetic contribution to A(q) results from 
exchange and short-range spin correlations, and is as- 
sumed to dominate a weaker, attractive electron-phonon 
contribution to the pairing interaction. The interaction 
depends on the correlation length, ^af, the incommen- 
surate wave vectors for the spin-excitation spectrum, 
Sx,y, and the electron-phonon and magnetic coupling 
strengths, Aop and Aaf. The quasiparticle dispersion 
relation is modelled by the tight-binding parametriza- 
tion for YBCO taken from Ref. The phonon- 

mediated interaction is assumed to be attractive in the 
spin-singlet, s-wave (Aig) pairing channel. Indirect evi- 
dence for an attractive sub-dominant s-wave interaction 
in the cuprates comes from ab-plane tunnelling mea- 
surements on Pr-doped YBCO which show strongly sup- 
pressed transitions, Tc < 20 K, and no ZBCP indica- 
tive of d-wave superconductivity jl3 |. For this pairing 
interaction and band structure we calculated the eigen- 
values Ar and eigenfunctions, nrf pf ) , of the linearized 
gap equation following Refs. [l^. The eigenfunc- 
tions determine the anisotropy of the pairing state in 
momentum space and form basis functions for the irre- 
ducible representations (T) of the point group, D4h- The 
eigenvalues determine the instability temperatures Tj, 
for pairing into states with symmetry dictated by the 
corresponding irreducible representation. The pairing in- 
teraction in Eq. ^ allows us to examine the sensitivity 
of the dominant and sub-dominant pairing channels to 
the parameters defining the model: ^af, Sx,y, and the 
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coupling parameters, Aep,af, and to consider the doping 
dependence via these parameters. We consider the one- 
dimensional, spin-singlet representations, Aig,A2g,Big, 
and B2g. The most attractive eigenvalues calculated for 
the two-channel interaction are shown in Fig. ^ as a 
function of 6. Note that all four channels have attractive 
eigenvalues for incommensurate spin-fluctuations i.e. for 
q vectors away from the (±7r/a, ±7r/a)-points. This is 
a robust feature of the model for ^af ^ 4 a. Secondly, 
the Big (d2;2_j,2) representation is the dominant pairing 
channel for 6 < 7r/2a. A third key point is that the 
leading sub-dominant pairing channels are nearly degen- 
erate over a wide range of values of the incommensurate 
wavevector, i5, for Aep = 0. Thus, weaker interactions like 
the electron-phonon coupling can play an important role 
in determining the sub-dominant pairing channel. For 
^af = 2a and aSx,y = f varying the strength of Aop allows 
us to tune the leading sub-dominant pairing channel; in- 
creasing Acp favors the Aig-channel, and it becomes more 
attractive than the B2g channel at Acp — O.lAaf. 
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FIG. 1: Eigenvalues and eigenfunctions for the singlet, ID 
irreducible representations of based on Eq. Q (a) Varia- 
tion of Ar with incommensurability parameter, S, for ^af = 2a 
and Aep = 0. (b) Basis functions for the dominant (Big) and 
the sub-dominant (B2g and Aig) pairing channels for aS = 
The black curve is the Fermi surface. 



The eigenfunctions, ?7r(p/), also depend on the param- 
eters of the pairing model: ?7Bigj VB2g and ?7A2g depend 
most sensitively on ^af and S^-^y, while r/Aig also depends 
on Acp. The calculated basis functions are shown in Fig. 

for Acp ~ 0; T^Big is well described by the simplest 
da;2_j,2 basis function, rid{(t>p) = a/2 cos (2 ^p), where 
is the angle p/ makes with respect to the [100] direction 
of the crystal. For the B2g channel the calculated basis 
function is not the simplest dxy basis function, ijd'{4'p) = 
-\/2 sin(2(/!)p), but a higher harmonic which to a good 
approximation is ?7B2g(0p) — 2 | sin(20p)| sin(6(/>p). The 
Aig basis function is highly anisotropic for Acp < O.lAaf 
but nearly isotropic for Acp > O.lAaf. The Aig basis 
function is to good approximation given by, '7Aig(0p) = 
(2/\/l — 2q! + Sa^) (a — (1 — a) sin(2(/)p) sin(6^p)) , where 
the parameter a varies between a = for Aep <C O.lAaf 
and a = 1 for Acp 3> O.lAaf. In general a varies with dop- 



ing and provides us with a convenient paramctrization of 
the effect of doping on the pairing eigenfunctions. 

The basis functions, 77r(p/), Tc, and the relative cou- 
pling strengths of the sub-dominant pairing channel are 
inputs to calculations of the surface phases and tun- 
nelling conductance for d-wave models of the cuprates. 
The surface OP, A(p/,R) = Ar(R)77r(p/), conden- 
sate momentum, p^, and excitation spectrum are ob- 
tained, selfconsistently, from solutions for the propaga- 
tor, spectral density and OP describing the supercon- 
ducting state. The pr opag ator obeys the quasiclassi- 
cal transport equations [iJl, and the OP, A(pj^,R), is 
determined by the anomalous quasiclassical propagator, 
5'(p/, R; e„), through the BCS gap equation, 

A(p/,R)=T^ [ rfV/A(p/-p});?(p^,R;e„). (2) 

The pairing interaction, A(p/ — p'y), is resolved in terms 
of its eigenfunctions and eigenvalues: A(p/ — Pj) = 
^P Arf?r(p/)?yf (p'y); and wc eliminate the interaction 
parameters, Ar, and cutoff frequency, Ec, in favor of the 
measurable instability temperatures, The for- 

mulation and application of the quasiclassical transport 
equations to surface states and tunnelling is described in 
Refs. IIEIEI 

The surface superconducting phase is sensitive to mag- 
netic fields that generate surface screening currents. 
These currents give rise to Doppler shifts of the low- 
energy surface excitations, eDoppio, ~ ■ Ps(R)j which 
modify the relative stability of surface phases described 
by different local symmetries, and may induce surface 
phases with broken T-symmetry. This type of Reld- 
induced transition leads to low-field nonlinearities in the 
field-dependence of the splitting of the ZBCP, S{H). The 
condensate momentum is given by the phase gradient and 
gauge- invariant coupling to the vector potential, p^ (R) = 
ICVi? — (e/c)A(R)). For magnetic fields applied along 
the crystal c-axis, Ps(R) is perpendicular both to c and 
to the surface normal n, and of order VfPs — A{H/Ho) 
where Ho = {c/e)A/vfX is the pairbreaking field of or- 
der of a few Tesla in the cuprates 0. The Fermi ve- 
locity, Vf and density of states, Nf, determine the zero- 
temperature penetration depth 1/A^ = (47re^/c^) A^/w| 
in the clean limit, and we take A/^o = 100 for YBCO. 

The conductance of a NIS tunnel junction measures 
an angle-average of the quasiparticle density of states 
(DOS) at the interface. In the limit of a low-transmission 
tunnel barrier (2? <C 1) the conductance is given to 
leading order in V by the standard tunnelling result, 
dl/dV = {1/Rn) J-^ d^pfV{pf)N{pf,Ks;eY) in the 
limit T <C Tc, where the integral is over trajectories with 
P/-n > 0, Rn is the interface resistance, iV(p/, R^; eV) = 
— ^Um ©■'^(p/, Rs; e = eV) is the angle-resolved local 
DOS of the superconductor evaluated at the surface, Rs 
and 0^ is the diagonal element of the retarded propa- 
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gator, 25^ ~ ©(icn e + 10+). The angle-resolved lo- 
cal DOS is folded with the trajectory dependent barrier 
transmission probability, I?(p^). The tunnelling barrier 
is modelled by a transmission probability that is maxi- 
mum for quasiparticles with Fermi momenta parallel to 
the surface normal (p/||n) and decreases with angle as 
V = I?o c"''^/'^''-' where (j) = cos~-'^(p/ • n) and (fic is a 
measure of active tunnelling trajectories. Qualitatively, 
thick junctions have narrow tunnelling cones, (j)c "C 7r/2, 
while thinner junctions correspond to a wider range of 
active tunnelling trajectories. 
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FIG. 2: The field dependence of the low-bias tunnelling con- 
ductance, dl/dV, calculated at different fields for a dx2_y2- 
wave superconductor in the presence of a sub-dominant 
Aig(a = 0.5) pairing interaction [panel a], and a sub- 
dominant B2g pairing interaction [panel b]. The pairing 
strength is too weak for a local T-breaking state to form spon- 
taneously at the surface. 



In Fig. 121 we show results for the tunnelling conduc- 
tances for the sub-dominant pairing channels with Aig 
and B2g symmetries. The surface DOS was calculated at 
fixed temperature just above a surface phase transition 
with broken T symmetry, T > . The conductances 
were calculated in the clean limit at T = O.I Tc. Note 
the nonlinear field evolution shown in Fig. |3i for the 
Aig(a = 0.5) channel. For the B2g channel the splitting 
of the ZBCP is suppressed at low fields. 

The field induced spfitting of the ZBCP, S{H), is lin- 
ear in H for pure d-wave pairing as shown in Fig. 
|31 The proximity to a broken T-symmetry phase of the 
form d -I- is enhances the splitting of the ZBCP and gen- 
erates a field dependence of the shift in the conductance 
peak that is nonlinear in H ioi T > Tg (Figs. |3i and 
^.). By contrast, the field splitting for the dxy chan- 
nel is suppressed at low fields; but onsets above a crit- 
ical field of order H*, = {T^^" /Ta)Ho (Fig. Et). The 
field dependence of the spfitting, 5{H), and its sensi- 
tivity to a sub-dominant order parameter, Ar, can be 
understood by considering the ABS spectrum for an 
OP that is constant in space, but breaks T-symmetry, 
e.g. A(p/) = Ad(p/) -f iAr(p/). The dx2_y2 compo- 
nent. Ad, changes sign along the scattered trajectory, 
while the sub-dominant component is unchanged. Thus, 



A(pj)* = — A(p/) on each trajectory. From the re- 
tarded propagator, (5^(p/,e), we obtain the ABS pole 
at an energy shifted away from the Fermi level, e;,(p/) = 
-sr Ar(p/) - V/ • p^, with sr = sgn[Ad(p/)Ar(p/)]. 




FIG. 3: Splitting of the ZBCP, 5{H), as a function of mag- 
netic field and strength of the sub-dominant pairing channel, 
for ij>c = 7r/4. Field-dependence of the sub-dominant order 
parameter at the surface. Sub-dominant channels: (a,b) Aig, 
(c,d) Bag. 

In the absence of the Dopplcr shift, a sub-dominant OP 
shifts the ABS. The magnitude and sign of the shift de- 
pend on the magnitude of the sub-dominant OP and the 
relative phase between Aj and Ar, sr = ±- The bound- 
state energies for trajectories, p/ and — p/, have oppo- 
site signs; this splitting for time-reversed states results 
in spontaneous surface currents carried by the bound 
states 0, . The two possible directions for the sponta- 
neous current reflect the two-fold degeneracy of the sur- 
face phase, i.e. Ad(p/) ± iAr(p/). In a magnetic field 
screening currents generate a Doppler shift of the sur- 
face bound states. For the sub-dominant s-wave OP the 
Doppler term may generate a shift of the same sign as the 
spontaneous shift, or the reverse depending on the polar- 
ity of the field and the relative phase of the subdominant 
OP. If the shifts are opposite in sign, the bound state en- 
ergies move toward the Fermi level with increasing field, 
until it becomes energetically favorable for the relative 
phase and bound state current to reverse sign. This re- 
versal occurs at a low field of order H* = {Ts/Tc)Ho for 
the sub-dominant s-wave channel as shown in Fig. ^jp. 
Also shown is the field-induced sub-dominant OP that 
onsets rapidly with field for T > Tg, and which is the 
source of the nonlinearity in d{H) shown in Fig. For 
T < Tg this splitting is asymmetric about H = 0, and 
hysteretic on a field scale of order H* . 

In the doping range where the leading sub-dominant 
channel is B2g, the field evolution is more complex. The 
angular dependence of the B2g OP, Ad^^y ~ sin 2((/)p — f ), 
implies that the bound state disperses through the Fermi 
level as a function of the trajectory angle, (j> ~ cos^^(p/ • 
n). For a [110] surface half of the incident trajectories 
correspond to ABS above the Fermi level and the other 
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half have their energies below the Fermi level. When 
added to their time-reversed partners, these states carry 
oppositely directed currents. Field-induction of the dxy 
OP for T > Ts produces a different behavior than pre- 
dicted for s-wave sub-dominant pairing; the magnitude of 
Ad^y remains small, and the corresponding bound state 
splittings are compensated by the Doppler term, until 
a threshold field is reached, again of order H*, beyond 
which the Doppler term cannot be compensated by the 
shift from the field-induced sub-dominant OP (see Fig. 
ISt-d). For T < Ts the subdominant OP produces split- 
ting of the ZBCP for H = 0. The low-field evolution 
is again different than that of the s-wave case; 5{H) 
decreases with H until a critical field of order H* is 
achieved. This field corresponds to the counter-moving 
branch of ABS dispersing through the Fermi level. Fur- 
ther increasing the field can drive the sub-dominant OP 
to very small values, in which case the splitting is given 
by the Doppler splitting and linear for H ^ Ho- 



a field- induced dx2_y2 -l-idxy surface state; there is a low- 
field threshold before a splitting appears, followed by an 
approximately linear increase in the splitting with H. 
At optimal doping the linear field dependence is recov- 
ered with no spontaneous splitting. Further increase in 
the doping level shows the nonlinear regime for slightly 
overdoped samples, which is accounted for by a attrac- 
tive s-wave subdominant pairing channel. The evolution 
with field and doping is systematically reproduced by a 
Big + iAig(a) surface state with anisotropic Aig compo- 
nent closest to optimal doping which evolves as the dop- 
ing is increased to a nearly isotropic Aig sub-dominant 
pairing state. The evolution of the ZBCP indicates that 
the pairing interaction changes with doping; the domi- 
nant pairing channels is Big, with a sub-dominant B2g 
component in the slightly underdoped regime which is 
overtaken by a sub-dominant Aig component in the over- 
doped regime. The cross-over occurs close to optimal 
doping. 
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FIG. 4: (a) Calculated field dependence of the conductance 
maximum as a function applied field for different values of 
doping parameter a. The model for the transmission prob- 
ability is = Vo expl- {4>/(j)cf], with (f)c = 12°. (b) 
Measured field dependence of the conductance maximum as 
a function of applied field for difi'erent values of doping as 
reported in Fig. 3 of Ref. 

Many of these features appear to be observed in YBCO 
tunnel junctions. We conclude with a comparison of the 
theoretical predictions of the surface ABS model for the 
tunnelling conductance, including the qualitative predic- 
tions of the two-channel pairing model, with the tun- 
nelling measurements reported in Ref. [lOj for the field 
dependence of the tunneUing conductance for YBCO. 
Calculated splittings of the ZBCP, S{H), are shown in 
Fig. and are compared with the experimental re- 
sults shown in Fig. Experimental observations of the 
field dependence of the conductance peak in underdoped 
films are in agreement with theoretical results based on 
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